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ABSTRACT 

The CARMA 1.3 mm polarization system consists of dual-polarization receivers that are sensitive to 
right- {R) and left-circular (L) polarization, and a spectral-line correlator that measures all four cross 
polarizations (RR, LL, LR, RL) on each of the 105 baselines connecting the 15 telescopes. Each receiver 
comprises a single feed horn, a waveguide circular polarizer, an orthomode transducer (OMT), two 
heterodyne mixers, and two low-noise amplifiers (LNAs), all mounted in a cryogenically cooled dewar. 
Here we review the basics of polarization observations, describe the construction and performance of 
key receiver components (circular polarizer, OMT, and mixers—but not the correlator), and discuss 
in detail the calibration of the system, particularly the calibration of the R-L phase offsets and the 
polarization leakage corrections. The absolute accuracy of polarization position angle measurements 
was checked by mapping the radial polarization pattern across the disk of Mars. Transferring the Mars 
calibration to the well known polarization calibrator 3C286, we find a polarization position angle of 
X = 39.2 ± 1° for 3C286 at 225 GHz, consistent with other observations at millimeter wavelengths. 
Finally, we consider what limitations in accuracy are expected due to the signal-to-noise ratio, dynamic 
range, and primary beam polarization. 

Keywords: instrumentation: interferometers; instrumentation: polarimeters; techniques: polarimetric; 
techniques: interferometric; polarization 


1. INTRODUCTION 


Millimeter-wave radiation from astronomical objects 
may be polarized by a number of processes, including 
many that involve magnetic fields (e.g., synchrotron ra¬ 
diation, Zeeman splitting, dust emission, etc.). This 
paper describes the design, calibration, and performance 
of the 1.3 mm (210-270 GHz) receiver system used for 
polarization measurements with the Combined Array for 
Research in Millimeter Astronomy (CARMA). CARMA 
is an 23-antenna aperture synthesis telescope; the 1.3 mm 
receivers are installed on the six 10 m-diameter and nine 
6m-diameter antennas. 


The CARMA system makes it possible to observe the 
polarized thermal emission from dust in protostellar cores 
with angular resolutions of 1-2", allowing one to study the 
magnetic field morphologies on scales of ^ 1000 AU in the 


nearest star-forming regions (Hull et al.||2013l 

Krumholz 

et al. 

2013 

IStephens et al.|20131|Hull et al.|2014 

IDavidson 

et al. 

2014 

IWright et al.|2014l|Stephens et al.|2014l|Segura- 

Cox et al.p015 ) or toward evolved stars (|Sabin et al.|2015 ). 


to measure Fara day rotation toward t he active galactic 
nucleus in 3C84 ( Plambeck et al.[2014 ), constraining the 


mode of accretion onto the black hole in this source, and 
in Very Long Baseline Interferometer experiments that 
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promise to probe the event horizons of nearby black holes 
( Johnson et al.|[^Q14 ). 

The plan of this paper is as follows. We first review the 
measurement of polarization via the Stokes parameters, 
and we motivate our choice of crossed circular feeds for 
the CARMA system. We then describe the design and 
performance of the key receiver components—broadband 
waveguide circular polarizers, orthomode transducers, and 
mixers—and discuss in detail the calibration of the R-L 
phase offsets and polarization leakage terms. We describe 
measurements of the absolute polarization position angle 
(PA) using Mars, and transfer this calibration to source 
3C286. Finally, we consider what limitations in accu¬ 
racy are expected due to the signal-to-noise ratio (SNR), 
dynamic range, and primary beam polarization. 


1.1. Stokes parameters 


It is conventional to characterize the polarization prop¬ 
erties of the incoming radiation field in terms of Stokes 
parameters. The Stokes parameters fully describe the 
characteristics of radiation that is fully polarized, par¬ 
tially polarized, or unpolarized. In an X-Y coordinate 
system where -\-x points North, -\-y points East,"^ and 
radiation propagates toward us along the -\-z axis, the 
Stokes parameters are given in terms of the time-averaged 
products of the complex voltages Ex and E y (Equations 
2.47a-2.47d in Rybicki & Lightman 1979[ Equations 1 


Hamaker fc Bregman||1996[ and Equations 4.19 ir 


Note that, as in Figure 4.8 of|Thompson et al.|(|2004|) the X 
and Y axes are rotated 90° counterclockwise relative to standard 
Cartesian axes. 
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Thompson et al.||2004|): 


I = 

(ExE^j + (EyEp) 

( 1 ) 

Q = 

{ExE*x) - (EyEp) 

(2) 

u = 

{ExEp) + (EpEy) 

( 3 ) 

V = 

-i{{ExEp)-{E*xEy)) , 

( 4 ) 


where denotes the complex conjugate of Ex- 

Alternatively, as described in more detail in the Ap¬ 
pendix, the Stokes parameters may be expressed in a 
circular polarization basis: 


/ = {ErEP) + {ElED 

( 5 ) 

Q = {ErED + {EpER 

( 6 ) 

U = -i {{ErED - (EpED) 

( 7 ) 

V = {ErEP) - {ErED . 

( 8 ) 


Note that Equations [T]-[4| and [SHSi ar e appropriate for 
a sing le dish telescope, ^s described in [Thompson et ah 
(2004), an interferometer like CARMA actually mea- 


sures Stokes visibilities from cross-correlations between 
antennas. For example, for telescopes m and n, Stokes 
Qmn = {ExmExn) “ {^YmEyn)• Thcsc visibilities are 
Fourier transformed to produce images of Stokes /, Q, U^ 
and V toward the astronomical source. These are maps 
of radio brightness in units of Jansky/beam. Stokes / is 
always positive, but Stokes Q, t/, and V can be negative 
or positive (but not complex). 

From the Stokes images, one can derive the fractional 
linear polarization Hi and position angle y: 


H/ 


+ f/2 


I 


1 U „ 

X = - arctan — , 0 < x < tt. 

1.2. Choice of crossed circular feeds. 


( 9 ) 

( 10 ) 


For polarization measurements with an aperture syn¬ 
thesis array, the receivers may be configured either with 
crossed linear {Ex and Ey) or crossed circular {Er and 
El) feeds. For observations at 1.3 mm wavelength, the sci¬ 
ence goals are mostly to measure weak linear polarization; 
few sources have appreciable circular polarization. Lin¬ 
ear polarization is derived from Stokes Q and U. With 
crossed linear feeds, Stokes Q (Equation is derived 
from the difference in power measured by the X and Y 
receivers. One must therefore take the difference of two 
large numbers {EyEy and EyEy) in order to measure 
a small number (Q), which demands that the X and Y 
gains be extremely stable and well calibrated. The gain 
stability requirement is relaxed considerably if crossed 
circular feeds are used, since then both Q and U are de¬ 
rived from cross-products of Er and Er- The magnitudes 
of these cross products are nearly zero if the source is 
weakly polarized, hence gain fluctuations are much less 
of a problem. For this reason, it is generally considered 
better to observe weak linear polarization with crossed 
circular feeds, and this is the approach that CARMA 
uses. For a more detailed discu ssion of the pr os and cons 
of linear and circular feds see Cotton (1998), as well as 
the summary (Section]^. 


1.3. Observing modes 

The 1.3 mm receiver system operates in three 
modes: single-polarization mode {LL or RR) dual¬ 
polarization {LL and RR simultaneously), and full-Stokes 
(l/L, RR, RR, RR). Depending on the correlator setup, all 
three modes can use up to the full 8 GHz bandwidth of the 
CARMA correlator (4 GHz per sideband), but with dif¬ 
ferent combinations of correlator bands and polarizations. 
See Table for the properties of the different observing 
modes in the wide-band (500 MHz/band) setup, which 
we used for observations of dust polarization. 

There is no benefit to using dual-polarization mode for 
wide-band (continuum) observing. Typically the RCP 
receivers have higher noise temperatures than the LCP 
receivers, so it is better to use the available correlator 
bandwidth to sample a wider frequency range in LCP, 
rather than a narrower range in both LCP and RCP. The 
dual-polarization mode offers higher sensitivity only for 
spectral line observations, where it is beneficial to get 
independent information from the second polarization, 
even if the receiver temperature for this polarization is a 
little higher. 


1 mm dual-polarization receivers 


SIS mixers 



Figure 1. The 1.3 mm dual-polarization receiver module, in¬ 
stalled inside the dewars on all the 10 m and 6 m antennas at 
CARMA. On the 10 m telescopes all these parts are cooled to 
4 K. On the 6 m telescopes, to reduce the heat load on the 4 K 
cryocooler stages, the two WBA13 amplifiers are connected 
to the 15 K cryocooler stages, while all other parts are cooled 
to 4K; the heat load from 15 K to 4K through the stainless 
steel SMA connectors is only 1-2 mW. 


2. HARDWARE 

Figure is a photo of the 1.3 mm dual-polarization re¬ 
ceiver module that is mounted in the cryogenically cooled 
dewar on each telescope. It includes a sin gle feed horn, a 
waveguide circular polarize r (Se ction 2.1), an orthomode 
transdu cer ( OMT; Section 2.2), two heterodyne mixers 
(Section |2.3| ), and two low-noise amplifiers (LNAs). The 
local oscillator (LO) and sky signals are combined using 
a mylar beamsplitter in front of the dewar window. 
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Table ^ CARMA 1.3 mm wide-band Observing Modes 


Mode 

Bands 

(#) 

BW/band 

(MHz) 

Pol. 

(#) 

SB 

(#) 

Total BW 
(GHz) 

Cross-corr. 

Channels 
(^ per band) 

Single-pol 

8 

500 

1 

2 

8 

RR or LL 

95 

Dual-pol 

4 

500 

2 

2 

8 

RR and LL 

95 

Full-Stokes 

4 

500 

2 

2 

8 

RR, LL, RL, LR 

47 


Table 1 

1.3 mm CARMA wide-band (500 MHz/band) observing modes. The total bandwidth of each mode is 8 GHz, and is equal to the 
number of correlator bands x the bandwidth per correlator band (BW/band) x the number of polarizations (Pol.) x the 
number of sidebands (SB). Observations with narrower bandwidths also are possible, and provide a greater number of channels. 


2.1. Waveguide Polarizer 


The CARMA dual-polarization receivers use waveguide 
polarizers, operating at temperatures of 4K, to convert 
incoming R and L polarized signals into X and Y lin¬ 
early polarized signals. These polarizers have lower loss 
and are more compact than quarter waveplates mounted 
outside the dewars. Since the polarizers are permanently 
installed inside the dewars, they must operate over the full 
210-270 GHz tuning range of the receivers. A polarizer 
consisting of a single quarter-wave retarder section would 
have much too narrow a bandwidth for this. However, it is 
possible to stack one or more half-wave retarder sections, 
rotated axially with respect to the quarter-wave retarder, 
to achieve broader bandwidth . Such multi-section po- 
larizers were first described by Pancharatnam (1955) in 
the context of birefringent wave plates. The action of 
the multi-section design is most easily visualized on a 
Poincare sphere (Figure]^. 

The CARMA polarizer is similar to the 2-section waveg¬ 
uide p olarizer developed by Kovac for t he DASI exper¬ 
iment (jLeitch et al.l 20021 Kovac||2004|). Whereas the 
(30 GHz) DAbl polarizer used dielectric vanes as retarder 
sections, the much higher frequency (230 GHz) CARMA 
polarizer uses sections of faceted circular waveguide as 
retarders. The guide wavelengths differ parallel and per¬ 
pendicular to the facets, causing a phase shift. The design 
of t he polarizer is described in det ail in CARMA Memo 
54 (Plambeck & Engargiola [2010 ). Figure shows the 
waveguide layout, and I'igurej^ gives the polarizer dimen¬ 
sions. The performance of the polarizer is characterized 
W the polarization leakages (see the discussion in Section 
1^. The leakages for ideal 1-, 2- and 3-section polarizers 
are shown in the left-hand panel of Figure the leakages 
expected for a 2-section polarizer when machining toler¬ 
ances are taken into account are shown in the right-hand 
panel. When one allows for machining tolerances, the 
performance of a 3-stage polarizer is no better than that 
of a 2-stage polarizer. 


Note that a multi-section polarizer cannot be flipped 
end to end—this is because rotations do not commute 
on the surface of a sphere. On the CARMA polarizers, a 
thicker flange with a dot marks the OMT end. 

As shown in Figure the polarizers are manufactured 
by electroplating gold and copper onto an aluminum 
mandrel. The copper electroform is machined to the 
correct diameter and soldered into a cylindrical hole in 
a brass shell. After machining is finished, the part is 


electroplated with gold and the aluminum mandrel is 
dissolved with sodium hydroxide to leave the finished 
waveguide. All electroplating steps were done by A.J. 
Tuck Co. (Brookfield, CT). 

We checked the performance of the polarizers by in¬ 
jecting a linearly polarized signal into the the horn end 
of each polarizer. By definition, this linearly polarized 
wave is the superposition of equal amplitude R- and L- 
polarized waves. Thus, for an ideal polarizer one would 
measure equal powers out the X- and T-ports of an OMT 
attached to the output port of the polarizer. Figure 
shows the results of these tests. The input signal was 
tuned from 205-270 GHz. Two sets of power measure¬ 
ments were made, with the OMT and two power meters 
rotated by 90° for the second test in order to average 
out differences in the OMT loss and the power meter 
calibrations. The theoretically expected power split for 
our polarizer design is shown in the lower left panel of 
the figure. 

The measured data for each polarizer were fitted to an 
analytic model (based on the one described in CARMA 
Memo 54) that allowed for errors in the angles and thick¬ 
nesses of the faceted waveguide sections. Smooth curves 
in each panel show the R and L powers predicted by the 
model for that polarizer. The model is a simplification in 
that it assumes that the facet depths of each section are 
uniform and equal for both retarder sections. However, 
while the fits suggest that the thicknesses of the faceted 
retarder sections normally were within the ± 0.0003" al¬ 
lowed tolerance, the angles of these retarder sections 
often were outside the ±0.2° tolerance. Although it was 
difficult to inspect the inside of the waveguide on the com¬ 
pleted polarizers, careful measurements with an optical 
comparator found angular offsets that were roughly con¬ 
sistent with the fitted results. The expected polarization 
leakages may be derived from the fitted dimensions; for 
most of the polarizers the predicted leakage amplitudes 
are less than 5% across the 210-270 GHz band. 


2.2. Orthomode transdueer 

After the waveguide polarizer converts incoming R and 
L polarized signals into X and Y in the circular waveg¬ 
uide, these linearly polarized signals are coupled into 
separate rectangular waveguides by an orthomode trans¬ 
ducer (OMT). As shown in Figure |8, the OMT employs 
a turnstile junction to split each m the two incoming 
polarizations into two opposite waveguide arms. These 
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Figure 2. Action of single- and 2-element circular polarizers visualized on a Poincare sphere. The principal axes of the sphere 
correspond to Stokes Q, and V. The goal is to transform R and L circularly polarized radiation, at the north and south 
poles, to X and Y linear polarization, at opposite ends of the Q-axis on the equator. Left: Passing RCP radiation through a 
single quarter-wave retarder with its principal axes aligned at 45° to the A-direction corresponds to a rotation on the sphere 
that places the central frequency (230 GHz, green) exactly on the equator, but spreads neighboring frequencies (210 GHz, red; 
260 GHz, blue) along an arc, leading to substantial ellipticity in their polarizations. Right: In a 2-stage polarizer the quarter 
wave section is followed by a half-wave retarder that rotates the center frequency back to the equator; the dispersion of this 
half-wave retarder nearly cancels chromatic errors produced by the quarter-wave section, thus broadening the bandwidth. 



Figure 3. Left: Waveguide layout. The polarizer converts an L-polarized signal entering the horn end to a T-polarized signal 
exiting the OMT end, while R is converted to X. Right: Outer view of the polarizer. A thick flange with a dot indicates the 
OMT end. The dot is aligned with the polarizer’s T-axis for serial numbers GP01-GP19, but (due to an unfortunate mixup) 
with the A-axis for GP20-GP36. 


arms recombine in T^-plane Y-junctions. 


blocks to join tightly at the tuning stub is difficult. 


The OMT is fabricated from four blocks that meet along 
a common edge. Initial tests of the design were made with 


a 20GHz scale m odel (Navarrini et al.|2005 Navarrini & 
|Plambeck||2006 ). To test the manutacturability at 1 mm, 
OMT's were then made by 4 independent machine shops 
(RAL, Custom Microwave, Protofab Inc., U of Arizona). 
These were tested by Alessandro Navarrini and Alberto 
Bolatto using a vector network analyzer at NRAO in 
Charlottesville, VA. De tails of these tests are reported 
in CARMA M emo 32 (Navarrini et a l.| 2QQ6 ), and are 
summarized in Navarrini & Flambeck| (|2(J(J6p. The test 
results showed excellent input return loss (better than 
-15 dB) and polarization isolation (better than -40 dB). 
The transmission losses were initially 1-2 dB, a little 
higher than expected, but improved by several tenths 
of a dB when gaps in the tuning stub at the base of 
the tunstile were filled with indium. Unfortunately, the 
tuning stub is split between the 4 blocks, and getting the 


Some of the prototype OMTs showed narrow resonances 
in their transmission curves. Simulations showed that 
these could be caused by small differences in the lengths 
of the opposite waveguide arms between the turnstile 
and the power combiner, perhap s due to misalignments 
of the 4 blocks (see Figure 19 in Navarrini et al.||2006). 
Energy coming from the turnstile junction reflects back 
from the power combiner if the two signals reaching it are 
not exactly 180° out of phase. Thus, standing waves can 
be set up between the turnstile and power combiner at 
frequencies where an integral number of half wavelengths 
fit within the waveguide arm. 


Subsequently, a set of 30 OMTs were machined by 
Protofab Inc. (Petaluma, CA). These brass parts were 
gold plated. All but one of the OMTs were tested on 
a demonstration N4252A network analyzer at Agilent, 
Inc. (Santa Rosa, CA), in 2010 January and August. 
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74.5°±0.2° 



0.5 


—0.0302±0.001 


-^1-IJ> 



R0.125C 



VIEW G-G 0.035±0.0003J 


Figure 4. Polarizer dimensions (in inches). 




200 220 240 260 


Figure 5. Left: Theoretical polarization leakage for 1-, 2-, and 3-section waveguide polarizers (red, green, blue curves) using 
faceted circular retarder sections. The dashed green curve shows the expected performance of a 2-section polarizer when 
machining tolerances are taken into account. Right: Simulated leakages for 200 polarizers with machining tolerances of ±0.0003^^ 
on the thickness of the faceted sections, and ±0.2° on the angles of the flats. The green curve shows the mean result that is 
reproduced (dashed) in the left-hand panel. 


We gratefully acknowledge Suren Singh for making these 
tests possible. Transmission measurements for the 15 
polarizers ultimately used on the telescopes are shown in 
Figure Average transmission losses are in the range 
0.4 to 0.7 dB. Narrow (^150 MHz FWHM [full width 
at half maximum]) resonances in the transmission, with 
losses of up to 2 dB, were often a problem, however. For 
the handful of OMTs measured both in January and 
August, the frequencies of the strongest resonances were 
identical. Often the resonances disappeared if we loosened 
the screws holding the 4 quadrants of the OMT together, 
and in some cases it was possible to retighten them in a 
sequence that maintained better performance. Loosening 
and tightening the screws changes the waveguide widths, 
hence the guide wavelengths, hence the electrical lengths 


between the turnstile junction and the Y-junction power 
combiners. The arms are about an inch long, so one 
expects resonances spaced by about 6 GHz, as observed. 

Occasionally an OMT that worked well at room temper¬ 
ature developed resonances when cooled to 4 K, probably 
because of differential thermal contraction between the 
screws and the 4 OMT quadrants. We were able to detect 
these resonances by measuring the receiver noise temper¬ 
ature across the 1-10 GHz IF passband over a range of 
LO frequencies. A lossy resonance in the OMT shows up 
as a narrow spike in the noise temperature, and this spike 
moves through the IF as the LO is stepped to different 
frequencies, forming what we termed “boat wakes” in the 
noise temperature plots. A good example is shown in 
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Soldered into 
waveguide flange 


Machined 


Copper 
electroplated 
onto mandrel 


P 

Aluminum 

mandrel 


Figure 6. Photo of the different stages of polarizer construc¬ 
tion. 


Figure Since it was not possible to adjust any of the 
screws when the OMT was at 4K, we chose simply to 
replace any OMT showing such behavior with a spare. 


2.3. Mixers 


The mixers use ALMA (Atacama Large Millimeter- 
submillimeter Array) Band 6 SIS (superconductor- 
insulator-superconductor) tunnel junctions fabricated at 
the University of Virginia by Arthur Lichtenberger. These 
are arrays of 4 SIS devices in series, offering greater dy¬ 
namic range and higher saturation power than a single 
SIS device. Most of our mixers were constructed using 
“C12,L56” or “C14,L56” devices from wafer 1489. 


Although at AL MA these devices are used in sideband- 
separating mixers (Kerr et al. 112013), the CARMA mixers 

DE 


are a simple double sideband (USB) design. An advantage 
of this design is that the local oscillator can be coupled 
into the mixer via a beamsplitter mounted outside the 
dewar. Although a DSB mixer folds together signals arriv¬ 
ing in the sidebands above and below the local oscillator 
frequency so that they appear at the same IF frequency, 
a phase-switching pattern applied to the local oscillators 
allows the correlator to separate the upper and lower 
sideband signals in cross correlation spectra. Unfortu¬ 
nately, this is possible only for signals that are common 
to a pair of telescopes, not for atmospheric noise that is 
independent for each telescope. Thus, DSB mixers allow 
each correlator segment to process twice its bandwidth in 
sky frequencies, but noise from both sidebands appears 
in the spectrum of each sideband. 


A magnet mounted behind the SIS junction is used 
to tune out Josephson tunneling through the insulating 
barrier, which otherwise adds noise and makes it difficult 
to tune the mixers stably. On ALMA an electromagnet is 
used for this purpose, so the field can be adjusted to cancel 
the Josephson tunneling very precisely; for simplicity, the 
CARMA mixers use just a tiny permanent magnet for 
this purpose. The magnet clings to a steel 4-40 set screw 
in the aluminum block attached to the mixer. The set 
screw may be adjusted (at room temperature) to position 
the magnet to get about the right field. 


One of the steps in optimizing the tuning of an SIS mixer 
is to adjust the LO power level. In the CARMA system 
the LO power cannot be adjusted independently for the 
RCP and LCP mixers. The LO is linearly polarized. It is 
coupled in through a beamsplitter. As shown in Figure 
a linearly polarized signal injected into the module wm 
be coupled about equally—but not perfectly equally— 
into the LCP and RCP mixer blocks. However, mixers 
sometimes prefer to operate at different power levels. It 
would be possible to install a rotatable quarter-wave plate 
between the LO waveguide horn and the beamsplitter 
to produce any elliptical polarization, which would offer 
the fiexibility of adjusting the power levels independently 
on the two mixers, but to keep the system as simple as 
possible this is not done. Instead, the procedure is to 
optimize the LO power for the LCP mixer only; the RCP 
mixer must live with whatever it gets. Thus the RCP 
noise temperatures at CARMA tend to be slightly higher 
than the LCP noise temperatures. 

Figure [TT] displays the DSB receiver temperatures as a 
function oFLO frequency for the mixers installed on the 
telescopes. These noise temperatures were measured in 
the lab using 295 K ambient and 77 K cold loads. They 
are based on total power measurements of the entire 1- 
9 GHz IF band using a broadband power meter, with no 
bandpass limiting filters. No corrections were made for 
optics losses in the dewar or for noise contributions from 
warm amplifiers. The best devices have DSB receiver 
temperatures of approximately 40 K over much of the 
210-270 GHz band. The median receiver temperature, 
shown by the thick black curve in the plot, is slightly 
greater than 50 K from 210-240 GHz. 

The average power and DSB noise temperature across 
the IF passband, at an LO frequency of 220 GHz, are 
shown in Figure The usable part of the IF band is 
between about 2 and 8.5 GHz. 


3. CALIBRATING POLARIZATION OBSERVATIONS 


One approach to calibrating polarization data is 
through the use of Mueller matrices. A Mueller matrix is 
a transfer function between the observed and the actual 
Stokes pa rameters. This is th e appro ach discussed, for ex - 
ample, by Sault et al. ( 1995|) and by Heiles et al. (2001b). 
We do not use this approacn because CARMA data typi¬ 
cally are analyzed with the MIRIAD^ software package, 
which breaks the calibration into separate passband, gain, 
R-L phase,^ and leakage steps. 


The passband and gain calibrations are handled by 
MIRIAD tasks MFCAL and SELFCAL. The passbands 
generally are stable for the duration of each 6-8 hour 
observation, and the gain phases of the R and L channels 
track each other very closely as well. The R vs. L 
gain amplitudes vary substantially only if one receiver 
is mistuned; generally amplitude variations are due to 
pointing or focus errors, which apply equally to the R 
and L channels. 


^ This font signifies the data reduction package MIRIAD itself; 
a task, p rocedure, or keyw ord within MIRIAD; or a Python task. 
MIRIAD ( [Sault et al.|1995| > is one of the standard data-reduction 
software packages tor (subjmillimeter-wave interferometry. 

® While we are actually measuring the R-L phases at CARMA, 
the quantity is known in MIRIAD as XYphase. 
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Figure 7. Test results for all the circular polarizers used on the telescopes. A linearly polarized signal (which can be represented 
as a superposition of LCP and RCP with equal amplitude) was injected into the “horn” end of the polarizer, and the normalized 
power emerging from the R (red) and L (blue) outputs of the OMT are plotted as a function of frequency. The power levels 
would be perfectly equal for an ideal polarizer; the expected power split for a 2-section polarizer of our design in shown in 
the lower left panel. The data for each polarizer were fitted to an analytical model that allowed for errors in the angles and 
thicknesses of the faceted sections. / is the depth of the facets (the target dimension is 0.0060''), 03 is the angle of the half-wave 
section relative to the Y axis of the OMT (the target value is 15.0°), and 02 is the angle between the quarter-wave and half-wave 
sections (the target value is 59.5°). 


The R-L phase calibration corrects for the phase dif¬ 
ference between the R and L channels on each telescope 
(i.e., the R-L phase) caused by delay differences in the 
receiver, underground cables, and correlator cabling. The 
R-L phase is not a single number, but is a function of 
frequency, due to fiber- and cable-length differences. That 
is, the same piece of the IF in different correlator sections 
can have completely different R-L phases. See Section]^ 
for a detailed discussion of R-L phase and position-angle 
calibration. 

The leakage corrections compensate for cross-coupling 
between the R and L channels, caused by imperfections 
in the polarizers and OMTs, reflections inside the dewars, 
and cross-coupling of IF signals in the analog electronics 
both preceding and within the correlator. See Section 
for a detailed discussion of leakage calibration. 


position angle x.^ The ordinary passband correction done 
with MFCAL analyzes only on the parallel-hand signals 
LL and RR; it does not solve for the R-L phase difference. 
In order to fit the R-L phase it is necessary to observe 
a linearly polarized source with a known polarization 
position angle. R-L phase calibration on astronomical 
sources is difficult at millimeter wavelengths, since most 
calibrators are weakly polarized and their polarization 
position angle varies on time scales of weeks or months. 
At centimeter wavelengths 3C286 is the usual polariza¬ 
tion calibrator. It is approximately 15% polar ized, and 
its position angle has been stable for decades (Perley & 
Butler||2013|). Unfortunately, its flux density at 23U GHz 
is only about 0.4 Jy. 


4.1. Grid calibration 


4. R-L PHASE CALIBRATION 

For circularly polarized feeds, the phase difference be¬ 
tween the R and L channels is what measures the position 
angle of an incoming linearly polarized signal. An R-L 
phase difference of 2% corresponds to a linearly polarized 


^ For example, a 180° rotation of the R-L phase leads to a 90° 
rotation of x- This is because of the circular-to-linear conversion: 
an R-L phase change of 90° will transform the radiation from 
linear (R and L in phase) to circular (90° out of phase); another 
90° change will transform the radiation from circular to linear (180° 
out of phase, with a position angle perpendicular to the original 
linear orientation). 
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Figure 8. Left: OMT design. A turnstile junction splits incoming signals of each polarization into two rectangular arms 
that recombine in a waveguide F^-plane junction. Right: The polarizer is constructed from 4 blocks that meet along a common edge. 


Not only is it impractical to calibrate the R-L phase 
using astronomical sources, but matters are further com¬ 
plicated because in the CARMA system there is not a 
single R-L phase difference for each telescope. Rather, 
there is a different R-L delay —hence, a phase slope as a 
function of frequency—for each correlator band. 

These band-dependent delays arise in a number of 
places in the signal path: (1) on the front end, the path 
lengths from the input port of the OMT to the R and 
L outputs differ by 0.1 inch.^ (2) There are chromatic 
effects in the polarizer and OMT that cause slight dif¬ 
ferences in the phase shifts fo r sig nals in the upper and 
lower sidebands—see Section [4^ for further discussion. 
(3) There are differences in the R and L fiber lengths from 
the receivers to the correlator. (4) In the block down- 
converters, there are cable-length differences between the 
upper (> 5 GHz) path (the path with the mixer) and the 
lower (<5GHz) path (without the mixer). 

To overcome these various difficulties, we use obser¬ 
vations of artificial linearly polarized noise sources to 
calibrate the R-L phase. The noise sources are created 
by inserting wire grid polarizers (Figure [I^ into the 
beams of the 10 m telescopes. With the grid in place, 
one linear polarization reaching the receiver originates 
from the sky, while the other originates from a room 
temperature load. Since the room temperature load is 
much hotter than the sky, the receiver sees thermal noise 
that is strongly polarized (Figure p^. 

By default, the master observing script performs grid 
observations approximately every 45 minutes during full 
polarization observing tracks. The calibration requires 
between 1 and 2 minutes: 45 seconds for the grids to move 
in and out of the beams, and 30 seconds of integration. 

® A distance of 0.1 in is a couple of wavelengths at 230 GHz. 
However, both the incoming radiation and the LO travel through 
this same length, and it is only the difference in phase shift at these 
two frequencies that matters. For an IF frequency of 5 GHz, this 
differential phase shift is ~ 20°. 


The grid data are identified by setting the purpose key¬ 
word to P in the data header. Grid observations are done 
at the end of phase calibration observations. Although 
the sky signal includes flux from the calibration source, 
this is negligible compared with the huge signal from the 
ambient load, and thus does not affect the measured R-L 
phase.^ 

Since the polarized noise is local to each telescope, it 
does not show up in cross-correlations with other tele¬ 
scopes. However, it leads to a strong LR autocorrelation 
(or “cross-auto”) signal, which is the cross-correlation 
of the L and R channels from a single antenna. The 
MIRIAD task XYAUTO averages together all the LR 
autocorrelation data for the grid observations in each 
dataset to create a channel-by-channel passband correc¬ 
tion for the 10 m telescopes. Rewriting the data with 
UVCAT or UVCAL applies these corrections. In the new 
dataset, the phases of all the A and L channels are equal 
on the 10 m telescopes. This means that a linearly polar¬ 
ized signal reaching the receiver with the same PA as the 
noise source will produce LR and RL correlations with 
phases of zero. One of these 10 m telescopes must then 
be used as the reference antenna for the regular passband 
correction performed with MFCAL. The passband correc¬ 
tion synchronizes the R phases on all telescopes with the 
R phase of the reference antenna, and the L phases on 
all telescopes with the L phase for the reference antenna. 
Since the R and L phases of the reference antenna are 
equal, then the R phases of all antennas equal the L 
phases on all antennas. 

Absolute position angle of the grids. MIRIAD 
computes the polarization position angle of a source as 

® For example, 3G279 has a brightness of ~ 10 Jy, and is about 
10% polarized. The antenna gains of the 10 m telescopes are about 
65 Jy/K, so in terms of brightness temperature, the flux from the 
3G279 is roughly 15 mK. This is two to three orders of magnitude 
smaller than the (polarized) temperature difference between sky 
and ambient load, which is of order 100 K even at low elevation in 
bad weather. 
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frequency (GHz) 

Figure 9. Transmission loss (S21) measurements of the OMTs used on the telescopes, made with an Agilent N4252A network 
analyzer. Losses in waveguide transitions to the OMT have been calibrated out. Red curves are the transmission from the 
circular input to the “short” arm of the OMT (rectangular waveguide port closer to the input); black curves are transmission to 
the “long” arm. The step size is 18.75 MHz. The lossy resonances are roughly 150 MHz wide at the half-power points. 


X = g' + evector, where q is the source’s parallactic angle^^ 
and the uu-variable evector is defined, for linear feeds, 
as the position angle of the A-feed relative to the local 
vertical. For circular feeds, evector is interpreted as the 
position angle x for which RCP and LCP radiation are 
in phase with one another, which is determined by the 
angle of the wire grid noise source. 

With the grids in place and the telescopes pointed at 
zenith, noise reaching the receivers is vertically polar¬ 
ized. However, because of the reflection off the tertiary 
mirror, a vertically polarized signal in the receiver cabin 
corresponds to a horizontally polarized (x = 90°) source 
on the sky (Figure [l^. Consequently, in datasets from 
CARMA the MIRIAD uu-variable evector is set to 90°. 

The parallactic angle is defined as the angle between two 
great circles, one through a celestial object and the zenith, the 
other (the “hour circle”) through the object and the celestial poles. 
The hour circle is analogous to a longitude line on a globe. In the 
zenith-object-celestial pole triangle, the parallactic angle is the 
angle at the celestial object’s corner. The parallactic angle is zero 
when the object crosses the meridian. 


4.2. Systematic limits to R-L phase accuracy 

The grid calibration of R-L phase is susceptible to a 
number of limitations that an astronomical calibration 
would not be. 

USB and LSB are averaged together. The autocor¬ 
relation averages together the R-L phase differences in the 
upper and lower sidebands. In normal cross-correlation 
spectra between antennas m and n, the LSB and USB are 
separated by demodulating the phase-switching pattern 
between the m and n local oscillators. This is not possible 
for the LR autocorrelation spectrum for a single antenna, 
since both the RCP and LCP channels use a common 
LO. 

The delays that contribute to the R-L phase difference 
are mostly at the IF frequency, and thus produce precisely 
the same phase shifts for signals in the LSB and USB. 
The only exception is the difference in delay from the 
beamsplitter to the RCP and LCP mixers. There is a 
0.1 in length difference through the OMT for the two 
polarizations, and the waveguide is slightly dispersive. 
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I.F. freq (GHz) 

Figure 10. Receiver noise temperature vs. IF frequency, mea¬ 
sured at 46 different LO frequencies, from 210 GHz (bottom) 
to 255 GHz (top) in 1 GHz steps. In this plot resonances in 
the OMT show up as sharp peaks that march across the IF 
band, forming “boat wakes.” In this example mixer 54-22 was 
mounted on the short arm of OMT 10. The handful of OMTs 
that showed this behavior were disqualified for use on the 
array. 



Figure 11. DSB receiver temperatures as a function of LO 
frequency for the mixers used on the array (except 4 that are 
not shown because their data were lost in a disk crash). These 
are based on total power measurements of the entire 1-9 GHz 
IF passband with a broadband power meter, with no band 
limiting filters. The median receiver temperature is shown by 
the thick black curve. 

so the delay is slightly different for the USB and LSB. 
A quick calculation shows that this effect is small: the 
phase shift of the LSB relative to the LO is at most 2° 


average CARMA receiver performance at 220 GHz (mean of 22 mixers) 



IF frequency (GHz) 

Figure 12. IF output power and DSB noise temperature as a 
function of IF frequency, from an average of 22 mixers. Noise 
temperatures were based on measurements of the IF power 
made with an Agilent E4407 spectrum analyzer, using a 3 MHz 
resolution bandwidth. 

greater than the phase shift of the USB relative to the 
LO.^^ The waveguide polarizer also is slightly chromatic, 
so the phase shifts through it will differ slightly for the 
upper and lower sidebands. 

Common noise. Noise that is common to both the 
R and L channels will generate an LR autocorrelation 
signal even when the wire grid is out of the beam. For 
example, noise radiated out the input of one mixer could 
be transmitted out through the OMT, polarizer, and 
feed horn, reflect back into the module from the dewar 
window (as the opposite circular polarization), and be 
coupled into the opposite mixer. If the R and L mixers 
had independent local oscillators with different phase 
switch patterns, these signals could be rejected in the LR 
autocorrelations; however the local oscillator is shared 
between the two polarizations in the CARMA receivers, 
so this is not possible. 

In order to derive the R-L phase with high accuracy, the 
signal from the polarized noise source must be much larger 
than the background level. To test this, we compared the 
LR amplitudes in typical 1 mm weather with and without 
the grids in place. The amplitudes are ^ 20 times higher 
with the grids in. This corresponds to an rms error of 
^1/20 radian, about 3°, which leads to an uncertainy of 
1.5° in position angle. Of course, in poor weather or at 
low elevation the uncertainty can be greater because the 
contrast between sky and ambient is less, so the polarized 
noise level is lower. 

Leakage. There is always some crosstalk between the R 
and L channels due to polarization leakage. As described 
in Section the leakages exhibit considerable frequency 
structure. Ideally one would solve for R-L phase and leak¬ 
age in an iterative way, but this is not easily accomplished 
with MIRIAD, since there is no simple way of applying 

The guide wavelength is Xg = Aq/^I — (Ao/2a)2, where Aq 
is the wavelength of the radiation in free space and a is the broad 
dimension of the waveguide. The dispersion is greatest at the low 
end of the band. Taking an LO frequency of 210 GHz and an IF of 
8 GHz, the phase shift is 36.2° for the LSB and 34.8° for the USB. 
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Figure 13. Left: Photo of a wire grid polarizer in the receiver cabin of a 10 m telescope, with the telescope pointed at zenith. 
The wires in the grid are vertical. Right: Photo of the tertiary mirror with the 10 m telescope pointed at zenith. Horizontally 
polarized radiation from the sky reflects off the tertiary such that it is vertically polarized at the receiver. The receiver with its 
wire grid rotates in elevation along with the telescope, so this correspondence is preserved at all elevations. 




Figure 14. A schematic showing how the wire grid polarizers 
produce a highly polarized noise source by reflecting ambient 
radiation into the receiver. 

leakage corrections on a channel-by-channel basis. This 
means that the leakages can introduce ripples into the 
R-L phase calibration. 

Phase-lock loop (PLL) jitter. Each correlator band 
uses separate digitizers for R and L. The digitizers are 
not run off of a single 1 GHz clock, but instead are run 
off of independent clocks that are phase-locked to a com¬ 
mon reference signal. One systematic effect that limits 
our ability to measure absolute position angle is rapid 
variation in R-L phase solutions caused by jitter (phase 
noise) in the 1 GHz PLL outputs. See Figure for plots 
of the R-L phase residuals (after correction), which show 
1-2° variations on very short (10 s) timescales. This effect 
should average out for most astronomical measurements. 

5. LEAKAGE CALIBRATION 

Leakage corrections compensate for cross-coupling be¬ 
tween the R and L channels, caused by imperfections in 
the receivers or crosstalk in the analog electronics that 



0 50 100 150 0 50 100 150 0 50 100 150 

Channels Channels Channels 


Figure 15. A sequence of cross-auto (LR) grid calibration 
data on antenna C3, taken at 10 second intervals, after the 
R-L phase correction has been applied. The 1-2° residuals 
are caused by phase jitter in the digitizer clocks. 

precede the correlator. Leakages are measured in terms 
of voltages (1% [or 0.01] leakage in voltage corresponds 
to 10~^ in power), and are defined in the following way 
( Thompson et al.||2004[ Equation 4.42): 

Vr = vrL- DrVl ( 11 ) 

^'l=^lLDlVr, ( 12 ) 

where v'^ and are the observed signals, vr and vr are 
the true signals. Dr is the leakage from L into R, and Dr 
is the leakage from R into L}‘^ Note that the leakages 
are complex numbers. 

Leakages are calibrated by observing a strong source 
(usually the gain calibrator) over a range of parallactic an¬ 
gle. The calibrator may be polarized or unpolarized. Since 
the telescopes have altitude-azimuth (alt-az) mounts, the 

In MIRIAD, Dx —>■ Dji^ and Dy —)■ Dr. 






























12 


Hull & Plambeck 


source’s polarization will appear to vary, in the frame of 
the receivers, as the telescopes track it across the sky, 
modulating the LR and RL cross correlation amplitudes 
in a predictable way. The component of the cross correla¬ 
tion amplitude that does not vary with parallactic angle 
must then be due to instrumental leakage. 

MIRIAD task GPCAL fits observations of the cali¬ 
brator in order to solve simultaneously for the antenna 
gains vs. time, the polarization leakages, and the source 
polarization. Given good weather and normal antenna 
performance, a 4-6 hr observation for which the source 
parallactic angle varies by more than about 60° yields a 
reliable leakage solution. 

Unfortunately, MIRIAD allows for only a single pair 
of leakage corrections and Dy for each telescope. To 
solve for channel-by-channel leakages requires one to run 
GPCAL many times, specifying each channel range in 
turn and saving the results. We have written a library 
of Python routines to handle this chore, and have ac¬ 
cumulated a library of leakage solutions covering many 
frequency ranges. 


Expected leakage amplitudes. Ideally, the polariza¬ 
tion leakages would be due exclusively to imperfections 


/eiy 

in the waveguide polarizers. Figure shows that we ex¬ 
pect leakage amplitudes of ^ 2-3% for polarizers that are 
within the expected dimensional tolerances. From the 
polarizer test data in Figure we inferred that many of 
the polarizers were not within these tolerances; however, 
the leakage vs. frequency expected from each polarizer 
can be computed if i ts dimensions are known, using the 
software described in Plambeck & Engargiola (2010). 


Figure shows the band-averaged leakages for all tele¬ 
scopes based on observations obtained from 2011-2015. 
In some cases dewars were swapped, so fewer datasets 
were used. The dashed curve on each plot shows the 
theoretical leakage expected for the particular circular 
polarizer on that telescope, based on the fits to the po¬ 
larizer dimensions shown in Figure In most cases the 
leakage amplitudes are substantially larger than expected 
from the polarizer alone, although for a few telescopes 
(Cl, C2, C8, Cl5) the polarizer curve appears to form a 
lower bound to the measured leakages. 


With higher frequency resolution, one finds that in most 
cases the leakages have substantial frequency structure, 
with periods of a few x 100 MHz to a few GHz. Figure pT| 
shows examples of this structure for C2, one of the best 
antennas, and for C13, one of the worst. The ripples in 
the leakages suggest that cross coupling of the unwanted 
polarization takes place via multiple paths with different 
delays. We discuss possible sources of cross-coupling in 
the sections below. 


Cross-coupling in the block downconverter. One 

source of leakage ripples is cross-coupling of IF signals 
in the correlator room, probably in the block downcon- 
verters. Each of the 8 analog downconverters assigned 
to a telescope obtains its input from a 4-way switch on 
that telescope’s block downconverter. The 4 inputs to 
the switch are the low-band RCP, low-band LCP, high- 
band RCP, and high-band LCP, where “low-band” is the 
1-5 GHz piece of the IF band, and “high-band” is the 5- 


9 GHz piece (which has been downconverted to 1-5GHz). 
Although nominally the switch provides 45 dB isolation 
between ports, this is not adequate if the RCP and LCP 
power levels differ substantially. For example, if the RCP 
power is 15 dB greater than the LCP power, then the 
relative level of RCP coupled into the LCP IF is -30 dB, 
which means that \Vr\ = 0.03 |Ul|- This RCP signal 
beats with the RCP voltage coupled via the polarizer, 
producing a 6% peak-peak ripple in Dr. Meanwhile, the 
relative level of LCP in the RCP IF is -60 dB, which 
produces just 0.2% ripple in Dr. 

Thus, cross-coupling in the block downconverter mani¬ 
fests itself as leakage ripples in one polarization, but not 
the other. An example of this behavior is shown in Figure 
[Tsl In this observation 2 correlator sections were centered 
at 2.5 GHz in the IF, while the other 2 were centered at 
7.5 GHz. Prominent leakage ripples are seen only in Dr, 
and only in the 1-5 GHz section of the IF band. Evidently 
the RCP and LCP power levels are more closely balanced 
above 5 GHz. 

Eigure shows that the ripples can be reduced by 
equalizing the RCP and LCP power levels. Eor these 
observations the RCP power into the block downconverter 
on telescope C9 was initially 9 dB greater than the LCP 
power. Installing a 10 dB attenuator on the RCP input 
to the block downconverter substantially reduced the 
ripple in Dr for the the low-band sections, and caused 
no deterioration in performance above 5 GHz or in Dr. 


Reflections in the receiver. The ripples that we at¬ 
tribute to cross-coupling in the block downconverters have 
periods of about 250 MHz. Many antennas, particularly 
the 6 m antennas (C7-C15), also have leakage ripples with 
a period of about 1 GHz (see the results for C13 plotted 
in Figure pT|), which corresponds to cross-coupling with 
a Insec delay. We hypothesize that this is caused by 
reflections inside the dewar. In the 6 m dewars the front 
of the feed horn is approximately 5 cm behind the dewar 
window, and the path lengths through the components 
are roughly 5 cm for the horn, 2 cm for the waveguide po¬ 
larizer, and 3 cm for the OMT. Thus, an RCP signal that 
reflected off the RCP mixer would travel back through 
a 15cm path to the dewar window (as RCP). Reflection 
off the window would convert it to LCP. It then would 
travel through another 15 cm path to the LCP mixer, 
accounting for a total delay of 30cm/c, or Insec. 


To test this hypothesis, we attempted to reduce the 
reflection from the RCP mixer in antenna C9 by changing 
the mixer bias (see Figure 20). Reducing the reflection 
off of one mixer should reduce the leakage ripples for the 
opposite polarization. For the data shown in red in Figure 
the RCP mixer was biased to 13 mV. This is above 
tEe superconducting energy gap, so no RCP astronomical 
signals would have been downconverted to the IF: and 
indeed, the RCP leakages (red curves, lower panel) are 
just noise. However, at this bias there should be a better 
impedance match between the SIS mixer and waveguide, 
reducing the magnitude of the signal reflecting back from 
the mixer. And, indeed, the LCP leakage ripples (red 
curves, upper panel) appear to be reduced, consistent 
with the reflection hypothesis. 


We then tried to reduce the reflection from the dewar 
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Figure 16. Leakage amplitudes vs. frequency for all telescopes, derived from 12 datasets obtained over a 3.5 year period. 
Each point shows the leakage amplitude for a single 0.5 GHz wide correlator window. Red points indicate Dr; blue points, Dl- 
The dashed curve in each plot shows the leakage expected from the circular polarizer on that antenna, based on the polarizer 
measurements shown in Figure Fewer data are plotted for some antennas because of receiver swaps. 


window (a double convex Teflon lens) by tilting it with 
a special clamp ring. It was physically possible to tilt 
the lens by only about 5°, and this had no apparent 
effect on the leakages. In another experiment, we at¬ 
tempted to worsen the reflection from the dewar window 
on C7—an antenna with relatively small leakage ripples— 
by installing a flat 0.010 in thick piece of Mylar just in 
front of the window. Again, however, there was no dis- 
cernable effect on the leakages. Finally, to confirm that 
the ripples originate in the dewar and not elsewhere in 
the telescope or correlator room, we physically swapped 
the dewars between antennas C9 and CIO. This swapped 
the leakage pattern between these antennas, confirming 
that the problem does originate in the dewar. 


Reproducibility of the leakage solutions. Since the 
CARMA receivers have no moving parts, the leakages 
are expected to be stable and reproducible. Figure 
presents a histogram of the scatter in leakage solutions 
from several different datasets obtained over a two-month 
interval. In this case the standard deviation of the Re 
and Im parts of the leakage terms was ^ 0.009; however, 
in data taken a few days apart the standard deviation 
can be as small as ^0.002. 

Factors limiting the reproducibility are reflections from 
the mixers—which will depend on their physical temper¬ 
ature and on their voltage and current bias—and the 
relative RCP and LCP power levels in the block downcon- 
verters, which will depend on the tuning and correlator 
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Figure 17. Leakage amplitudes vs. frequency for antennas C2 and C13 over the frequency range 214-224 GHz. The frequency 
resolution is typically 50 MHz. In most cases there is excellent reproducibility in the leakages obtained months or years apart. 
Dashed curves show the leakage amplitude expected from the circular polarizers on these telescopes. 










Figure 18. Leakage amplitudes for antennas C2 and C3 plotted as a function of sky frequency, derived from a dataset where two 
correlator sections were positioned above 5 GHz in the IF (outer two columns, “LSB-high” and “USB-high”), and two below 5 GHz 
(middle two columns, “LSB-low” and “USB-low”); note that the lowest sky frequency corresponds to the highest IF frequency in 
the LSB. Black curves are Dr and blue curves are Dl- Large ripples are evident only in Dl, and only below IF frequencies of 
5 GHz (the “low” bands), indicative of cross-coupling in the block downconverters. 


setup. 

There is an ambiguity in the leakage terms because 
they always appear in pairs in expression s for the LR and 
RL cross-correlations (Sault et al.||l995 For example, 
the observed RL* cross correlation bet Wen antennas m 


The ambiguity does not occur for VLBI observations where 
sources are observed at different parallactic angles from different 
observatories. 


and n is given by: 

i'^Rm'^Ln) ~ ( T DRjyiV{vRji T ) 

(13) 

= i'^Rm'^Ln) ^Rmi'^Lm'^Ln) ^Rm"^Rn) 

(14) 

= i'^Rm'^Ln) ^Ln) 

This expression is unchanged if one adds an arbitrary 
complex number c = a + to all the Dr terms, and 
subracts its conjugate c* = a — jh from all the Dr terms: 


{DRm^a^jb)L{DLn-a^jby =DRm^Dl^. (16) 
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frequency (GHz) 

Figure 19. Dl leakage amplitude vs. frequency for antenna C9 before (black) and after (red) attenuating the RCP 
input to the block downconverter by 10 dB. The frequency setup is similar to that in Figure Lowering the RCP 

power level in the downconverter reduced coupling of RCP into LCP, hence reduced the magnitudes of the leakage rip¬ 
ples m. Dl for the two “low-band” correlator sections. The Dr leakages were unchanged, as were the leakages on all other telescopes. 



Figure 20. Leakage amplitudes as a function of frequency 
for antenna C9 derived from 9 different datasets. For the data 
shown in red the voltage bias on the RCP mixer was set above 
the SIS junction’s superconducting energy gap so that the 
junction behaved as a resistive load. This reduced reflections 
of RCP signals from this mixer, which lowered the magnitude 
of the LCP leakage ripples. The RCP leakages for this dataset 
were then just noise, since there was no signal coming from 
the RCP receiver. 


Figure 21. A histogram of the differences of 4 leakage so¬ 
lutions from their mean. The values encompass the real and 
imaginary components for Dr (leakage from LCP into RCP) 
and Dl (leakage from RCP into LCP), for 6 correlator win¬ 
dows and all 15 telescopes. The leakage solutions are from 02 
September 2012, 25 September 2012, 25 October 2012, and 
30 October 2012. The standard deviation is 0.0089. 


By convention, MIRIAD chooses this complex off¬ 
set c such that '^'Re{Djirn) = and 

"^lm{Djirn) = As a consequence, the 

absolute leakages may change if some antennas are miss¬ 
ing from the data set. When plotted in the complex 
plane, Dr and Dr have a reflection symmetry across the 
imaginary axis for most antennas (see Figure [2^. That is, 
Re(Di^) = —Re{DL) and Im(Di^) = Im(DL). Thus, the 
MIRIAD convention causes the average of the imaginary 
parts of the leakages to be close to 0. If an antenna like 
C3, for which the leakages have a large imaginary offset, is 
omitted from the solution, then the imaginary parts of all 
the other antennas will increase by about 0.01. One must 
therefore be cautious when comparing leakages derived 
from different data sets. 

6. SYSTEMATIC LIMITATIONS TO POLARIZATION 
MEASUREMENTS 


There are several effects that limit the accuracy of po¬ 
larization measurements made with CARMA. (1) The 
accuracy is limited by the signal-to-noise ratio; noise bias 
leads to an overestimate of the polarization fraction in 
the low-SNR case. (2) There is a dynamic range limit: 
for a bright source that is weakly polarized, errors in the 
leakage calibrations can give false detections of polariza¬ 
tion. (3) The absolute position angle accuracy depends 
on the grid calibrations. (4) For extended (non-point-like) 
sources, the polarization varies across the beam. These 
effects are discussed in the sections below. 

6.1. Signal-to-noise limitations: debiasing polarimetric 

images 

Polarization measurements have a positive bias because 
the polarization P = ^JQ‘^ is always positive, even 

though the Stokes parameters Q and U from which P is 
derived can be either positive or negative. This bias has 
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Figure 22. Leakages plotted in the complex plane. Dr (red) and Dl (blue) for four 500-MHz wide bands centered at 224.5, 
225.0, 225.5, and 226.0GHz, from data set cl217.2D_2303c279.iniriad.2. Typically Dr and Dl have mirror symmetry about 
the imaginary axis. 


a significant effect in low-SNR measurements, i.e., when 
P < 3crp, where ap is the rms noise in the polarization 
maps. (The rms noise values in the Q and U maps are 
generally comparable, such that we set ap ~ gq ~ crp). 
The bias can be taken into account by ca kulating the 


bias-c orrected polarized intensity Pc (e.g., Vaillancourt| 


a 

SNR 


|2QQ6| see also [Naghizadeh-Khouei fc Clarke||1993| tor 
discussion of the statistics of position angles in low SN 
measurements). 

The probability density function (PDF) for the observed 


by the Rice distribution ( 

Killeen et al.| 1986, 

Bl, and Vaillancourt||2006 

Equation 6): 


PDF(Pe|P,ap) = ^/o 


exp [-(p2 + P‘^)/2al] . 

dp J 

(18) 

Thus, given the observed polarization P, one can calculate 
the true polarization Pc by finding the maximum (i.e., 


the most probable value) of the PDF in Equation 18 


For very significant polarization dete ctions (P > 5(7p) 


the simple high-SNR limit is valid (see |Vaillancourt||2006 
Equation 12): 




+ [/2 


(19) 


PDF(P|Pe,ap) = T/o exp [-(p2 + p2)/24] . 

(17) 

However, to calculate the debiased intensity of the true 
polarization Pc one needs the opposite: PDE(Pc|P, ap). 
If one assumes a uniform prior for the true polarization 
Pc, then by Bayes’s theorem PDE(Pc|P, ap) is the same 
as Equation pT} 


However, for low-SNR detections (P < 5ap), things are 
not so simple, and one must use Equation [Ts] to calculate 
the debiased polarization intensity. 


The position angle x and uncertainty dx (calculated us¬ 
ing standard error propagation) of the incoming radiation 
are 


X 


1 

= - arctan 
2 



( 20 ) 


A = 


1 ap 

2 P ■ 


( 21 ) 
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I.e., for a detection Pc = 2crp, the uncertainty in the 
position angle is 1/4 of a radian, or ± 14°. 

6.2. Limitations due to leakage uncertainties. 

The accuracy of polarization measurements also is lim¬ 
ited by uncertainties in the polarization leakage correc¬ 
tions. These uncertainties affect bright sources as well as 
weak ones; thus they provide a dynamic range limit for 
polarization measurements with CARMA. 

To test the effects of leakage variations on the solutions 
of polarization position angle and fraction, we performed 
simulations using UVGEN. As discussed above in Sec¬ 
tion the Re and Im parts of the leakages vary by up 
to a tew X 0.01 from track to track. Thus, we set up the 
simulations to probe the effects of Gaussian random er¬ 
rors in the leakages with standard deviations between 0.01 
and 0.03. The simulated uv data included thermal noise 
consistent with typical 1 mm weather at CARMA: system 
temperatures of 300 K and opacity r 0.3. The data 
comprised 8 GHz of continuum bandwidth, and included a 
point source with a ffux of 1 Jy and at an elevation of 30°. 
We used actual leakage terms from CARMA polarization 
data as the initial set of values. 

We performed two tests, examining how variations in 
leakage terms affect the PA and the polarization fraction 
(%) of (1) a snapshot observation (for example, when 
measuring the polarization of the passband calibrator, 
which is only observed for 10 min), and (2) a longer, ^ 6 hr 
observation of the calibrator or the science target. 

The simulations proceeded as follows: 

• Used UVGEN to simulate leakage-free data: either 
a 10 min snapshot, or a long track covering ±3 hr 
around transit 

• Copied “real” leakages (from a typical 1 mm 
CARMA polarization track) into the data file 

• Rewrote the data to apply the leakages, thus cor¬ 
rupting the data 

• Copied a negated version of the above leakages into 
the dataset after varying the leakages’ Re and Im 
parts using Gaussian random errors with standard 
deviations between 0.01 and 0.03 

• Used UVFLUX to find the % and PA of the source 

• Ran a Monte-Carlo test for point sources with 0%, 
1%, and 10% polarization, and calculated the stan¬ 
dard deviation of the polarization fraction and PA 
for each value of the leakage error 

Copying negated leakages into the dataset corrects the 
corrupted data exactly. However, copying a set of leak¬ 
ages whose Re and Im parts were modified slightly result 
in slight differences in the calculated PA and %. Per¬ 
forming these simulations gave us a handle on how robust 
our calculated PA and % values were to typical leakage 
variations. 

See Figur e!^ for the results. Unsurprisingly, the scatter 
in PA and % is larger for snapshot observations than 


for full tracks. Additionally, the scatter in PA is much 
larger for weakly polarized sources. Looking at the plots 
of the scatter in %, we conclude that the instrumental 
polarization caused by the typical leakage errors of ^0.01 
is ^0.3% for a snapshot observation and ^0.1% for a 
longer track. The corresponding variations in position 
angle are ^8° and ^3° for a 1% polarized source. 


Note that very small variations in the leakages can cause 
drastic changes to the measured polarization position 
angles if a source is very weakly polarized {Pc/1 ^ 0.5%). 
For example, when analyzing CA RMA data toward the 
protoplanetary nebula CRL 618 (see Sabin et al.|2014) the 
position angles varied by up to 90"^ from night to night 
depending on how exactly we reduced the data. The 
scatter in the Re and Im parts of the leakage solutions 
on the two nights was ^0.01, which is standard; however, 
in this case these slight differences led to significantly 
different Q and U ma ps. This was not a pro blem for the 
TADPOL sources (see|Hull et al.||2013l|2014|), which were 
on average at least a tew percent polarized, and which 
tended to be much fainter sources where the ability to 
detect polarization was limited by the SNR instead of by 
dynamic range. 


We therefore urge caution when interpreting the po¬ 
sition angles of sources with polarization fractions of 
< 0.5%. 


6.3. Absolute position-angle accuracy 


Even for a strongly polarized source, there is an uncer¬ 
tainty in PA due to the absolute accuracy of the R-L 
phase calibration. As discussed in Section]^ polarized 
noise sources in the cabins of the 10 m telescopes are used 
to measure the R-L phase. One of the 10 m telescopes 
(usually Cl, since it has relatively small leakages) is used 
as the reference antenna for the passband calibration. 
This transfers its R-L phase calibration to all the other 
telescopes, including the 6 m telescopes, which are not 
equipped with polarized noise sources. 


Observations of an astronomical source with a stable, 
well defined PA are required in order to check the accuracy 
of the PA measurements. At centimeter wavelengths, 
3C286 is the usual polarization calibrator; its position 
angle has been stable for decades (jPerley & Butler||2013). 
However, the PA of 3C286 increases slowly with frequency, 
from X = 33° at A > 3.7 cm to x = 36° at A 0.7 cm. At 
millimeter wavelengths, Agudo et al. (2012) measured 
X = 37.3 ±0.8° at A 3 mm and x = 33.1 ±5.7'' at A 1.3 mm. 
ALMA commissioning results at A 1.3 mm gave x = 39° 
(Stuartt Corder, priv. comm., 2013). 


Fortunately, observations of the polarized thermal emis¬ 
sion from a rocky planet or satellite provide an absolute 
standard by which the PA of 3C286 can me measured. 
The planet must be well resolved by the synthesized beam. 
The polarization is expected to be radial with respect to 
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Figure 23. Errors in PA and fractional polarization as a function of leakage errors. 


the center of the planet’s disk.^"^ 


The radial polarization pattern has be en seen before in 
for example, observations of the Moon ([H eijnsfe ^Drake 


19631 Davies fc Gardner 1966, White fc (Jogdell 1973p ; 


Heiles &: Drake (|l963|) found that the maximum polariza¬ 
tion of the Moon was ^ 20% when observed at A 21 cm. 
The Moon is much too large to fit into the primary beams 
of millimeter wave telescopes like CARMA or ALMA, 
however. Mars is a better choice for these telescopes. 


Figure shows 225 GHz polarization observations of 
Mars made with CARMA on 04 May 2014. The planet 
was 14.2" in diameter and the synthesized beam was 
2.9" X 1.7". Observations of Mars, 3C279, and 3C286 
were interleaved over a 3 hour long track. Leakages 
were derived from the 3C279 data. A careful analysis of 
the polarization vectors (Figure [2^ shows that they are 


The radiation that we receive is a mixture of thermal emission 
from the planet that is transmitted through its surface, and 3K 
background radiation that is reflected off this surface. If the surface 
is tilted with respect to our line of sight, as it is near the limb of the 
planet, the transmission and reflection coefficients are functions of 
the radiation’s polarization direction. The transmission coefficient 
is greater, and the reflection coefficient smaller, for radiation that 
is polarized parallel to the plane of incidence, which is the plane 
defined by the incident and reflected rays. (In fact, at Brewster’s 
angle the reflection coefficient for this polarization drops to zero). 
Thus, emission from the planet is preferentially polarized parallel 
to this plane, i.e., radially. 


skewed by about 1.8° from the radial direction. 

PA depends on the choice of reference antenna. 

Unfortunately, a different choice of reference antenna 
yields a slightly different absolute position angle. As 
shown in column 5 of Table for the 04 May 2014 data 
the offset ranges from -0.56° for C5 up to 8.7° for C6. 
These discrepancies are unlikely to originate from leakages; 
C3 has the worst leakages (magnitudes of order 10%), but 
its PA deviation is no larger than average. The deviations 
may originate in the optics between the wire grids and 
the primary mirrors of the 10 m telescopes, but further 
tests should be made. 

For each choice of reference antenna we also computed 
the PA of 3C286. Column 10 of Table |2] shows the 3C286 
position angles after subtracting the corresponding Mars 
radial deviation. The answers are surprisingly consistent, 
with an average PA value of 39.9°. Using another leakage 
solution derived from 3C279 observations two days earlier, 
on 02 May 2014, gave a similar result: 40.0° for the PA 
of 3C286. 

Table summarizes Mars polarization observations 
on two different days: 04 May 2014 (see Table ^ and 
17 Jan 2014. The difference in the final corrected PA 
for 3C286 between the two days gives some indication 
of the uncertainty in the calibration. For the January 
observation Mars was observed over a narrow range in 
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Figure 24. Mars polarization map generated using Cl as the 
passband (and hence R-L phase) reference antenna, based 
on data set cl217.2D_2303c279.miriad.2. Mars is 14.2'' in 
diameter and the synthesized beam (blue ellipse) is 2.9" x 
1.7". The Stokes I intensity is shown by black contours 
at 15, 20, 25, 30, and 40Jybm“^. The peak Stokes I flux 
density is 47 Jy bm“^ and the maximum polarized intensity is 
1.4 Jybm“^. Line segments indicate polarization orientation 
from the CARMA data; segment lengths are proportional 
to polarization intensity. The polarization of sphe rical, solid 
bodies like Mars should be radial (see Section 6.3 and Figure 


25). 


parallactic angle, and the planet’s disk was less well- 
resolved by the synthesized beam, so it is tempting to 
assume that those data are less reliable. O n th e other 
hand, primary beam polarization (Section 6.4) might 
have skewed the results of the May observations, when 
the planet’s disk was larger. Taking the average of the two 
measurements, we conclude that the polarization position 
angle of 3C286 at 225 GHz is approximately 39.2° ± 1°. 


USB — LSB position angles. For Faraday rotation 
measurements, one measures the difference in a source’s 
polarization position angle in the LSB and USB. Slight 
differences in the R-L phase difference are expected at 
these two frequencies because of chromatic effects in the 
waveguide polarizers. It is not possible to derive the USB 
and LSB R-L phases from the wire grid noise source, 
since the sidebands cannot be separated in noise that is 
local to each antenna. 

Column 4 in Table |3] shows the difference in the Mars 
radial deviation for the USB and LSB. Depending on 
the choice of antenna used for the R-L phase calibration, 
this value ranges from 0.12 to 2.81°. Column 8 shows 
the difference in PA measured on 3C286 for the USB 
and LSB. For reasons that are still unclear, there is a 
systematic difference of about 1.2° between the Mars and 
3C286 offsets. We have higher confidence in the 3C286 
values, since 3C286 is a point source. 


6.4. Primary-heam Polarization 


If a source is not perfectly on-axis, then additional 
distortions are introduced across the primary beams of 
the telescopes. To check for these variations in the in¬ 
strumental polarization, we observed BL Lac (a bright, 
highly polarized quasar) at 16 offset positions, eight of 
which were 12" and eight of which were 24" from the field 
center.Our aim was to characterize the “beam squint” 
and “beam squash.” 

Beam squint arises when the feeds of the telescope are 
tilted with respect to the optical axis of the primary 
reflector, causing the LCP and RCP responses to be 
slightly offset from the symmetry axis. Beam squint 
is discussed in Chu fc Turrin (1973); Adatia & Rudge 
(|1975|);|Rudge & Adatia| (|l978|). Squint manifests itself as 
a double-lobed pattern in Stokes V, which is the difference 
of RCP and LCP (see Equation [^. 

Beam squash is caused by differences in the beam widths 
of the orthogonal polarizations. This is caused largely by 
the slight difference in reflectivity of the two polarizations 
off the surface of the curved parabolic reflector. Squash 
manifests itself as a four-lobed “cloverleaf” pattern in 
the linearly polarized St okes Q and Lf maps. This phe¬ 
nomenon is discussed 


Napier ([l994[ [199^; th e term 


beam squash” was coined in Heiles et al.| ( |20Ula ). 


Maps of both squint and squash have been made by 
Rob ishaw (|2QQ8 ) for the Green Bank Telescope (GBT) 
and Heiles et al. (2001a) for Arecibo. 


Additionally, there is a large body of work on wide- 
held polarimetric calibration and imagi ng with the Karl 
G. Ja nsky Very Large Array (VLA). |Uson & Cotton 
( 2008) primarily discuss Stokes V (squint) calibration 
Cotton & Perley (2010) discuss both off-axis circular 
(squint) and linear instrumental polarization, and present 
a model for correcting for both. They measure the off-axis 
instrumental polarization for various frequencies using 
the method we employ at CARMA, where all antennas 
move to various offset positions simultaneously. 

A more robust method of removing polarization arti¬ 
facts is beam holography, where one or more “reference 
antennas” remain pointed at a strong source while the 
remaining antennas move in a raster pattern, thus mea¬ 
suring the primary beam o f each dish in all four Stokes 
parameters IQU V (see, e.g., Corder & Wright 2Q06[ Lamb 


& Corder|2008 who map the total power [btokes 7 J beam 
at CARMA). After mapping all four beams for all anten¬ 
nas, one can then use those beam models to correct the 
full-Stokes data. 


Beam squint. While the CARMA system is almost 
always used for measuring linear polarization, we are still 
able to make maps of Stokes U, and are thus able to mea¬ 
sure squint. Beam s quint is normally charac terized using 
the squi nt angle Rudge fc Adatia||1978 Equation 25; 
see also Napierp994 ): 


The analysis presented below only uses the data taken 12" 
from the field center. 
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Figure 25. Analysis of deviation from radial of the Mars polarization data shown in Figure The red circle indicates the 
diameter of Mars. Radial deviations are computed only at pixels where the polarized intensity is > 0.3 Jybm“^. The color scale 
on the deviation map runs from -10 to 10°. Values are weighted by polarized intensity when computing the mean deviation. 




. A sin Oq 
= arcsin —— 
47rF 


A^o 

47rF ’ 


( 22 ) 

(23) 


where A is the observing wavelength, Oq is the angular 
offset between the feed and the telescope’s optical axis, 
and F is the focal length of the primary reflector. The 
approximation in Equation is good when the angular 
offset Oq is small, which is usually the case. 


The CARMA 6 m antennas have a Cassegrain design, 
with feed horn o ffset by Aaz = 7.126^ and Ael = 1.91' 

with an effective 


Plambeck 


2000). Using Equation 23 


ocal lengtn of F' = 1186 inches, the theoretical squint 
produced by these offsets should be incredibly small: 
^ 0.002". The measured squint is significantly larger 
0.45"; see Eigure 26, right panel). The feeds in the 
10 m antennas are on-axis (and should thus be squint- 
free); however, there are multiple reflections in the optical 
path, including three curved mirrors with ^ 90° bends 
that are all in different planes. The odd number of off-axis 
reflections may be the cause of the 10 m squint, which is 


even worse than the 6 m squint 2.16"; see Eigure 26 


left panel). 

In order to fit for the squint we solved for Stokes V at 
each (12") offset position and for each 30 s integration 
on BL Lac. Note that we applied the gain and leakage 
corrections solved at the center position only. We cal¬ 
culated the offsets in the frame of the dish (in azimuth 
and elevation) by de-rotating the offset position (in RA 
and DEC) by the parallactic angle x associated with the 
given integration. In order to derive the squint, we fit the 
data to the difference of two circular Gaussians G 2 — Gi \ 

Gi = exp(— [(az — -h (el — Ael/2)^] /(2cr^)) (24) 

G 2 — exp(— [(az + + (el + Ael/2)^] /(2cr^)) , (25) 


where a = EWHM/(2y2Mp), and EWHM is the full- 
width at half-maximum of the primary beam at A 1.3 mm 
(30" for the 10m antennas and 56" for the 6 m antennas); 
Aaz and Ael are the vertical and horizontal components 
of the offset between the RCP and LCP beams, and 
the angle of the beam offset in the frame of the dish 
d = arctan (Ael/Aaz) — 90°, measured east from north 
(or counterclockwise from vertical, in the reference frame 
of the antenna). We assume that the amplitudes of the 
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(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

R-L cal 


Mars radial dev (deg) 



3C286 PA (deg) 


3C286 PA (deg) 

antenna 

LSB 

USB 

USB-LSB 

DSB 

LSB 

USB 

USB-LSB 

DSB 

corrected 

Cl 

0.69 

3.50 

2.81 

1.77 

40.71 

42.25 

1.54 

41.62 

39.85 

C2 

4.96 

6.20 

1.24 

5.29 

44.92 

45.07 

0.15 

45.13 

39.84 

C3 

4.21 

4.86 

-0.53 

3.97 

44.15 

43.51 

-0.64 

43.95 

39.98 

C4 

5.73 

6.87 

1.14 

5.74 

45.63 

45.51 

-0.12 

45.69 

39.95 

C5 

-0.14 

-0.02 

0.12 

-0.56 

39.71 

38.62 

-1.09 

39.29 

39.85 

C6 

8.23 

9.74 

1.51 

8.70 

48.19 

48.80 

0.61 

48.63 

39.93 


Table 2 

Deviation from radial of the Mars polarization orientations, and 3C286 position angles, as a function of the R-L phase 
calibration antenna. Lower sideband (LSB; 217.75 GHz), upper sideband (USB; 232.25 GHz), and double sideband (DSB; 
225 GHz) values are given. For reasons that aren’t clear, the DSB values are not necessarily the mean of the LSB and USB 
values. The last column shows the 3G286 DSB position angles corrected by the Mars DSB radial deviations. Leakages were 
derived from data taken on 04 May 2014, during an observation with wide parallactic-angle coverage. 


Date 

Mars diameter 

Beam size 

Mars PA dev, DSB (deg) 

3C286 PA (deg) 




Cl C2 C3 C4 C5 C6 

corrected 

17 Jan 2014 

bo 

2.49 X 2.15" 

0.54 4.96 2.54 4.39 -1.83 4.54 

38.5 

04 May 2014 

14.2" 

2.92 X 1.71" 

1.77 5.29 3.97 5.74 -0.56 8.70 

39.9 


Table 3 

Average PA deviations in measurements of Mars, and 3C286 position angles. Same as Table but for two separate observations 
on 17 Jan 2014 and 04 May 2014. The difference in the final corrected PA for 3G286 is unexplained, but may be because the 17 
Jan track was a short, snapshot observation and the 04 May track had wide parallactic-angle coverage. 


Gaussians are identical, and that the two beams are offset 
by equal and opposite amounts (±Aaz/2 and ±Ael/2) 
from the pointing center. 

Beam squash. Telescopes with native circular feeds like 
CARMA should see squash in both Stokes Q and U, which 
are both combinations of the cross-polarizations RL and 
LR (see Equations and Q . For telescopes like the GBT 
and Arecibo with native linear feeds, it should be easier 
to see squash in Stokes U, which is calculated using the 
cross-polarizations XY and YX (see Equation]^ and is 
thus unaffected by AX and YY gain variations tnat can 
plague Stokes Q (see Equation [^. Maps of be am squash 
have been made in Stokes U for the GBT dRo bisha^ 
2QQ8|) and i n both Stokes Q and U for Arecibo (|Heiles| 
et al.|2(J01a[ note that for Arecibo the gain variations had 
a minimal effect, allowing excellent maps of both Stokes 
parameters). 

We see no evidence in the Q or U CARMA maps for 
a quadrupolar squash pattern, which should have twice 
the frequency of the squint as one moves around the 
dish. However, variations in Q and U do lead to squint¬ 
like variations (i.e., a single crest and trough) in the 
polarization fraction and position angle. 

Variation in polarization fraction and position an¬ 
gle across the primary beam. While we saw no evi¬ 
dence of squash in the data, we did see squint-like behavior 
in Q, U^ and the PA and polarization fraction derived 
from them. This behavior is not expected theoretically, 
but is seen at Arecibo and discusse d in various publica- 


tions and technical memos including 

Heiles (|1999); 

Heiles 

et al. (2001a 2003): 

; Heiles & Crutcl 

ierl2D05f. 


We fit a sinusoid to find the peak-to-peak variation in 
the calculated PA and polarization fraction of both the 
10 m and 6 m dishes. The 6 m dishes showed very little 
variation in PA and polarization fraction, whereas the 
10 m dishes showed more (possibly because observations 
were further out into the primary beam). The 10m 
antennas had peak-to-peak variations of ^ 8° in PA and 
^ 0.02 in polarization fraction (0.02 represents a ^ 25% 
variation in the mean of the BL Lac polarization fraction). 
See Figure!^ 

7. SUMMARY 

In this overview of the 1.3 mm dual polarization receiver 
system built for CARMA, we described the design and per¬ 
formance of the key hardware components—circular polar¬ 
izers, orthomode transducers, and mixers—and discussed 
the calibration of the system for polarization observations, 
focusing particularly on the calibration of the R-L phase 
offsets and the polarization leakage terms. The leakages 
were found to exhibit considerable frequency structure, 
which we attribute to reflections inside the dewars and to 
cross-coupling in the analog block downconverter, which 
precedes the correlator. 

We discussed limitations on the accuracy of polarization 
measurements made with this system, including the ef¬ 
fects of signal-to-noise, leakage uncertainties, and primary 
beam polarization. The absolute accuracy of polarization 
position angle measurements was checked by mapping 
the radial polarization pattern across the disk of Mars. 
Transferring the Mars calibration to 3C826, we measured 
a polarization position angle for 3C286 at 225 GHz of 
39.2° ± 1°. 
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Figure 26. 6 m and 10 m beam squi nt, fit u sing data positions offset by 12'' from the pointing center. The solid curve is the the 

best-fit squint model (see Equations 24 and 25). 






Figure 27. Variations of PA and polarization fraction across the 6 m and 10 m primary beams, fit using data from positions 
offset by 12" from the pointing center. The solid curves are sinusoids fit to the data to find the peak-to-peak variations. 


Lessons learned. In retrospect, we are not sure that 
the choice of circularly polarized feeds was a good one. 
In a dual polarization system with circularly polarized 
feeds, signals that reflect off a mixer can bounce back 
from the dewar window with the opposite polarization, 
introducing frequency-dependent cross-coupling between 


the two polarization channels. Accurately characterizing 
these ripples in the polarization leakages is extremely 
time-consuming, since it requires observations of a strong 
calibrator over a wide parallactic angle range. Although 
the leakage ripples in the CARMA system tended to 
average out over the full passband, and were not a severe 












TADPOL Data Release 


23 


limitation for continuum polarization observations of dust 
or extragalactic synchrotron sources, these ripples did 
hamper our ability to search for polarization in spectral 
lines. 

In contrast, in a dual polarization system with lin¬ 
early polarized feeds,reflections would not cause cross¬ 
coupling between the two polarizations {XY and YX), 
but instead would create ripples in the receiver bandpass 
{XX and YY)^ which are easily calibrated with a brief 
observation of a bright continuum source. Of course, as 
discussed in Section crossed-linear feeds have their own 
problems. Specifically, in order to detect weak linear 
polarization with crossed linear feeds, the receiver gains 
must be extremely stable. The newest interferometer with 
crossed-linear receivers is ALMA. The promising contin¬ 
uum polarization results obtained thus far with ALMA 
suggest that gain stability is not a significant limitation; 
further tests of the ALMA system will reveal whether 
there are other issues unique to crossed-linear systems 
that were not encountered at CARMA. 
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APPENDIX: CONVERTING FROM LINEAR TO CIRCULAR 
STOKES PARAMETERS 

Here we will work through the details of converting the 
Stokes parameters from their commonly seen linear forms 
to their circular forms. As we described in Section in 
an X-Y coordinate system where -\-x points North, -\-y 

We note that, if desired, a beamsplitter can still be used to 
couple the local oscillator into the mixers in such a system, simply 
by rotating the plane of polarization of the LO with a waveguide 
twist or half wave plate so that it is oriented at 45° to the OMT 
axes. 


points East, and radiation propagates toward us along the 
-\-z axis, the Stokes parameters for crossed-linear feeds 
are: 


I = {ExEy + {EyE*y) 

(26) 

Q = {ExE*x) - {EyE*y) 

(27) 

U = (ExE^) + {E*xEy) 

(28) 

V = -i{{ExE*Y) - {E*xEy)) , 

(29) 


where denotes the complex conjugate of Ex- These 
equations are consistent with any number of textbooks 


and publications, includ ing Equations 2.47a-2.47d in|Ry- 
bicki fc Lightman||1979[ Equ ations 1 in|Hamaker V Breg- 


man||1996[ Equations 4.19 in |Thompson et al.||2QQ4[ and 


others. 


In radio astronomy Er is defined using the IEEE (Insti¬ 
tute of Electrical and Electronics Engineers) convention: 
the phase (j) of the sinusoid Ey lags the sinusoid Ex by 
90°. Thus, a right-circularly polarized wave (one with 
pure Er) can be represented in the following way: 


Ey = Ex^ * 
= -lEx 


(30) 

(31) 


As defined and discussed in Contopoulos fc Jappel (1974); 


(1996|);|IEEE|(|1997|), this lag results 
irized (RCF) radiation whose phasor 


Hamaker & Bregman 
in right-circularly polarized (RCF) radiation whose phasor 
rotates counterclockwise as viewed by the receiver. 


Additionally, in Contopoulos & Jappel (1974) the lAU 
deemed Stokes V to be positive if the signal has net 
RCP. To check this, we plug Ex = and a lagging 

Ey = into the crossed-linear equation for Stokes 

V (Equation , and we find that indeed, Stokes V is 
positive for RuTP radiation. 

To convert from linear to circular we use the linear-to- 
circular coordinate tr ansform reported in Equation 17 of 
H amaker et al.|(|1996|), and furt her elucidated in Section 3 
m Hamaker V Bregman ( 199 6|) (note that we have chosen 
the upper sign convention in [Hamaker fc Bregmau||1996 
which is the convention used in IHamaker et al.lll996l as 
well as in the references they quote): 


Ca = 


V2 



(32) 


Hamaker & Bregman (1996) explicitly refer to the above 
matrix as circular-rl, clearly defining the order of their 
matrix (top row ^ R, bottom row ^ L). We can then 
convert from [Ex^Ey] to [Er^Er]: 


' Er 

1 

■ 1 i' 

'Ex' 

.El_ 

“Vi 

l-i 

Ey 


(33) 


Er and Er are thus: 
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Er 


El 


-^{Ex + iEv) 
^ {Ex - iEx) 


(34) 

(35) 


As a check, we substitute Ey = —iEx (Equations 30 
and 31, which describe how Ey lags Ex in RCP radiation) 
into me above expression for El. We find that El = 0, 
as expected. 


One can also invert matrix to obtain the matrix 
required to convert from circular back to linear: 


Ex 

1 

1 l' 

~ Er 

Ey 

"71 

—i i 

.El _ 


this yields expressions for Ex and Ey: 


Ex 


Ey 


-^{Er + El) 
{-iEr + iEL) 


(36) 


(37) 

(38) 


We can substitute these expressions for Ex and Ey 
(Equation s and i nto the crossed-linear Stokes pa¬ 
rameters (EquationspO}]^ . We find the crossed-circular 
Stokes parameters toEe: 


I = {ErE*L + {ElED 

(39) 

Q = {ErED + {E*^El) 

(40) 

U = -i {{ErED - {E*rEl)) 

(41) 

V = {ErED - {ElED ■ 

(42) 


Note that occasionally the Stokes parameters are 
fined as ^ X the expressions listed in Equations 


39 


de- 


42 


(e.g., Equations 1-4 in Roberts et aL||1994 ). When tms 


the case, one can express the tour cross-correlations of R 
and L in terms of the Stokes parameters: 


ErE}, = I + V 

(43) 

ElEI =I-V 

(44) 

ErEI =Q + iU 

(45) 

E}iEl =Q-iU . 

(46) 
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